
CLAIMS 



What is claimed is: 



1. 



A method for determining a geomagnetic transmission function for the transmission of a 



population of particles, said particles having one or more rigidities, to a satellite in a known earth 
orbit, under known geomagnetic conditions, comprising the steps: 

(A) dividing said orbit into a plurality of sequential steps comprising at least a first step and 
a last step, each sequential pair of said steps being connected at a point, each of said points 
corresponding to a position of said satellite along said orbit, and specifying a value for Kp and a 
value for Dst to specify said geomagnetic conditions; 

(B) for a plurality of said points, selecting a plurality of arrival directions for the arrival of 
said particles at said satellites, wherein each of said arrival directions represents either an allowed 
or a forbidden trajectory of one of said particles; and 

(C) for each of said arrival directions at each of said plurality of said points, determining 
whether said arrival direction represents an allowed or a forbidden trajectory by tracing a path of said 
particle to arrive at said point with said arrival direction, until said particle path intersects a boundary 
of earth's magnetosphere, thereby indicating an allowed trajectory, or until said particle path 
intersects either earth's atmosphere or earth's surface, thereby indicating a forbidden trajectory, 
wherein said tracing of said path is performed by integrating in the time domain the Lorentz equation 

December 31, 1997 (1:59PM) 1 



F = m/d^/ ( _ t/) = _Q vxB 



wherein F is the force vector acting on said particle, m is mass, y is 1/(1 - vVc 2 )" 1/2 , v is said particle's 
velocity vector, t f is the travel time of said particle, Q is the charge of said particle, and B is the 
magnetic field vector acting on said particle, wherein B is given by 

B = B ^(r, t') + B tsyg(Kp, r, t') 

wherein B lGRF is the International Geomagnetic Reference Field promulgated by the International 
Association of Geomagnetism and Aeronomy, and B XSYG is the modified Tsyganenko field given by 
the sum of B (T) + B (T) + B (T) + B (RC) + B (RC) + B ( RC > + B < c ) + R ( c ) + J\ (Q + 

ine bum Ol U Xmg T r> Ym g ^ °Zmg ^ ^Xmg ^ ^Ymg . ^ °Zmg ^ ^Xsm ^ ^Ysm ^ tt Zsm + 

B Xsm CM) + B Ysm (M) + B Zsm (M) , wherein coordinates in the solar magnetospheric system are denoted sm, 
and wherein coordinates in the solar magnetic coordinate system are denoted mg; and wherein 

B Y m = Q,yz r ; 



b,co - ^ 
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c, + c h 
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wherein W(x,y)= 0.5 



x- x„ 



, 1_ [(x-x 0 ) 2 + Df ; 



xli + -/ D;J 




wherein Q T = 



W(x,y) 



£ r S T 



Si<4" 3,-p*f X-j- S-j 



SLl 



wherein z, = z-z s (x,y, iff), 



5 (x,y,(/) = 0.5 tan ^(x+ R c - A /(x+ R c ) 2 + 16 
-Gsin^y 4 (y 4 +L 4 y )~ 1 , 



wherein S T RC = -J/7 2 + (a TRC + ^ TRC ) , 

£t,rc = -y/z r + D T RC , 
D T = D 0 + 8y 2 + r AO) + / A(x)» 



wherein C, = -98.72 when Kp = 0, 0 + , C, = -35.64 when Kp = V,l,l\ C, = -77.45 
when Kp = 2\2,2\ C,= -70.12 when Kp = 3",3,3 + , C, = - 162.5 when Kp = 4,4,4 + , 



December 31, 1997 <1:59PH) 



3 



C, = -128.4 when Kp * 5", 

wherein C 2 = -10014 when Kp = 0, 0 + , C 2 = -12800 when Kp = 1\U + , C 2 = 
-14588 when Kp = 2",2,2 + , C 2 = -16125 when Kp = 3 ,3,3 + , C 2 = -15806 when Kp 
= 4~,4,4 + , C 2 = -16184 when Kp z 5", 

wherein a T = 13.55 when Kp = 0, 0 + ,a T = 13.81 when Kp = r,l,l + , a T = 15.08 when 
Kp = 2\2X, a T = 1 5.63 when Kp = 3",3,3 + , a T = 1 6. 1 1 when Kp = 4",4,4\ a T = 1 5.85 
when Kp £ 5", 

wherein D 0 = 2.08 when Kp = 0, 0 + , D 0 = 1 .664 when Kp = 1\1 ,1 + , D 0 = 1 .541 when 
Kp = 2-,2,2 + , D 0 = 0.9351 when Kp = 3",3,3 + , D 0 = 0.7677 when Kp = 4",4,4 + , D 0 = 
0.3325 when Kp * 5", 

wherein 1^ = 9.084 when Kp = 0, 0 + , R c = 9.238 when Kp = V,\,l + , R, = 9.609, 
when Kp = 2",2,2 + , R== 8.573 when Kp = 3",3,3 + , R,. = 10.06 when Kp = 4',4,4 + , 
= 10.47 when Kp ;> 5", 
wherein L y = 1 0 R E , 



Qrc xz p 




wherein Q RC = ^>C 5 l KC A S K cK^Ac) 
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D RC = Do + YrcMx)+Y .hi(x) 

h TRC = 0.5[l + x(x 2 + L 2 TRC )- 1/2 ], 

h, = 0.5{1- (x+ 16)[(x+ 16) 2 + 36]~ 1/2 }, 



C 5 (Dst) = - 1 0220+408.5 -Dst 

B XYZ (C) = C 3 (F\ y<z + F^ +C 4 (F + X)yiZ -F- x>yz ) 5 wherein 



p± I -c±rc± 



W c (x,y) 



x 1— k 



S^S^Cz+R,)] y 



P* W c (x,y) 
r, = — — ; — + 



x-^r- L + y — - 



1 



V <?x dy ) S^CzIRt) 



S^KziR.^ + xVy 2 ]^, 



W c (x,y)=0.5 



1- 



x- x 



0c 



[(x-x 0c ) 2 + L 2 xc ]^ 



x(i + y 2 / D2 yc)"'; 



Bx (M) = e ^[c 6Z cos i|r + (C 7 + C g y 2 + CjZ 2 ) sin 

B Y (M) = e^C^yz cos i|r + (C„y+C 12 y 3 +C 13 yz 2 ) sin and 

Bz (M) = eX /Ax [(C|4 + Cisy 2 + C]6Z 2 ) cos ^ + (Ci?z+ CigZy 2 + C]9Z 3) sin 

wherein C 6 through C 19 are given by: 
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= 0,0 + 


= r,i,r 


= 2',2,2 + 


= 3",3,3 + 


= 4",4,4 + 


* 5" 


^6 


1 CI 1 
1 .5 1 J 


z.i lo 


1 A/1 1 

2.041 


1 1 O 1 

J.lol 


1 /CAT 

3.00/ 


A AAA 

4.oyo 


r 1 


i\ 1 a 
51.10 


1^ &A 
55.04 


/lO A£L 

42.40 


An ^a 
4/. 50 


^ 1 1 A 

51.10 


/I A A A 
49.09 




-0.0/404 


A A*7/| 1 

-U.U /41 


A ATA1 1 

-0.0/01 1 


-0.132/ 


A 1 AAA 

-0.1006 


a An 1 

-0.0231 


L 9 


-.0 / /04 


A 1 AO 1 

-O.lOol 


A 1 Clft 

-0.15 fy 


A 1 OA/1 

-0.1 o04 


A 1 An 

-0.1927 


A 1 CA 

-0.1359 


p 


0.U05305 


A AAOOO/1 

0.003924 


A AA/1 A*7Q 

0.0040/© 


0.013o2 


A AO T C5 

0.03353 


A A 1 AO A 

0.01989 


c n 


1 i in 

-1.129 


1 /1C1 

-1,451 


1 1 A 1 

-1.391 


1 /lOO 

-1.455 


1 ^> AO 

-1.392 


O O AO 

-2.298 


^12 


0.001663 


A AAOA1 

0.00202 


A AA1 C5 

0.00153 


0.002962 


A AA 1 C A A 

0.001594 


A A A A A 1 1 

0.0049H 


/~i 

C, 3 


A AAAAOO 

o.ooo9oo 


A AA 111 
0.00111 


A AAATT7 

0.000727 


A AAAOA1 

0.000897 


A A AO A 1 A 

0.002439 


0.00342 1 


Cm 


ion 


21.37 


21.50 


H 1 A 

22.1 r 4 


oo /i 1 
22.41 


O 1 *7A 

21.79 




A A1A1 O 

-0.03015 


A A /I C CI 

-0.04567 


A A /I 1 AA 

-0.04199 


A A /I f\f\C 

-0.04095 


A A A AO r - 

-0.04925 


A i\C A A 1 

-0.05447 


^16 


a nooofi 

-o.Oiozy 


A AC1 Ol 

-0.05382 


-0.06523 


-0.09223 


A 1 1 O 

-0.1 153 


A 1 1 A A 

-0.H49 




-0.1283 


-0.1457 


-0.6412 


-1.059 


-1.399 


-0.2214 


^18 


-0.001973 


-0.002742 


-0.000948 


-0.001766 


0.000716 


-0.01355 


C 19 


0.000717 


0.001244 


0.0002276 


0.003034 


0.002696 


0.00H85 



and wherein L y =l0.0, D x =l3.0, L RC =5.0, L T =6.30, Yt = 4.0, 6=0.010, y { =l.O, R T =30.0 5 

Xo c =4.0, L xc 2 =50.0, and D yc =20.0, 
and thereby determining whether said particle's trajectory intersects either the boundary of earth's 
magnetosphere, thereby indicating an allowed trajectory, or intersecting earth's surface or earth's 
atmosphere, thereby indicating a forbidden trajectory. 



2. The method of claim 1, wherein said plurality of arrival directions comprises one or more 



randomly or pseudorandomly selected arrival directions. 
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3. The method of claim 1, further comprising the step: 

(D) for each particle rigidity, determine what fraction of particles of that rigidity will be 
transmitted to said satellite. 

4. The method of claim 1 , wherein said plurality of points comprises a complete set of points 
for said orbit. 

5. The method of claim 1 , wherein said plurality of points comprises a set of points for a portion 

f • 

of said orbit. 

6 . A method for determining, for a given particle environment outside of earth' s magnetosphere, 
what portion of a population of particles having one or more rigidities making up said particle 
environment will be transmitted to a satellite in a known earth orbit, comprising the steps: 

(A) performing steps (A) through (C) o f clainO , thereby computing a geomagnetic 
transmission function for said population of particles; and 

(B) applying said geomagnetic transmission function to said population of particles, thereby 
determining what portion of that population of particles will be transmitted to said satellite. 

7. A method for determining, for a given particle environment outside of earth' s magnetosphere, 
what portion of a population of particles having one or more rigidities making up said particle 
environment will be transmitted to a satellite in earth orbit, comprising the steps: 

(A) prompting a user to specify an earth orbit; 
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(B) determining whether said orbit is among a group of preselected orbits, each of said orbits 
in said preselected group of orbits having an associated predetermined geomagnetic transmission 
function for a range of particle rigidities, each of said predetermined geomagnetic transmission 
functions having been prepared in accordance with claim A ; and 

(C) for the case wherein said orbit is among said preselected group of orbits, applying said 
predetermined geomagnetic transmission function for said orbit to said particle environment outside 
earth's magnetosphere. 

8. The method of daimj^wherein said preselected group of orbits comprises a quiet shuttle 
orbit having a 450 km altitude and a 28.5° inclination, a disturbed shuttle orbit having a 450 km 
altitude and a 28.5° inclination, a quiet space station orbit having a 450 km altitude and a 51.6° 
inclination, and a disturbed space station orbit having a 450 km altitude and a 51.6° inclination, 
wherein said quiet shuttle orbit has a geomagnetic transmission function given by FIG. 4, wherein 
said disturbed shuttle orbit has a geomagnetic transmission function given by FIG. 5, wherein said 
quiet space station orbit has a geomagnetic transmission function given by FIG. 6, and wherein said 
disturbed space station orbit has a geomagnetic transmission function given by FIG. 7. 

9. The method o f claim j^ further comprising the step: 

(D) for the case wherein said orbit is not among said preselected group of orbits, performing 
a step selected from the group consisting of (a) returning an error message to a user, and (b) 
computing a geomagnetic transmission function for said orbit, in accordance with the method of 
claim 1. 
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10. A method for determining a flux of a species of solar ions having a specified atomic number 
between 3 and 92, and a specified kinetic energy, for a satellite in a near earth orbit, comprising the 
steps: 

(A) specifying an atomic number for solar ions for evaluation; 

(B) specifying a kinetic energy for said solar ions; 

(C) specifying a baseline model for said flux of solar ions, wherein said baseline model is 
selected from the group consisting of a worst day model, a worst week model, and a peak flux 
model; 

(D) in the case wherein said specified atomic number is greater than 20, selecting iron as an 
elemental spectrum model; 

(E) in the case wherein said specified atomic number is less than or equal to 20, selecting 
oxygen as an elemental spectrum model; 



(F) looking up a value for an elemental breakpoint, wherein said elemental breakpoint is a 
function of said elemental spectrum model and said baseline model, and wherein said elemental 
breakpoint is selected from the table: 



^ — Baseline model 
Elemental ^^^^^ 
spectrum model 


worst day or peak flux 


worst week 


iron 


24.23 MeV/nuc 


19.90 MeV/nuc 


oxygen 


15.94 MeV/nuc 


12.89 MeV/nuc 



(G) in the case wherein said baseline model is said worst week model, and said elemental 
spectrum model is iron, and said kinetic energy is greater than 127.93 MeV/nuc, calculating an 

unsealed flux, wherein said unsealed flux equals A sp * ( En/ Me % uc ) " Gsp , wherein A sp = 3 . 1 68 1 4x 1 0 6 
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(cm 2 sr MeV/nuc)' 1 , (En/ Me /4, c ) is said kinetic energy, normalized to be dimensionless, and G sp 



= 2.861; 

(H) in the case wherein the condition recited in step (G) is not satisfied, and wherein said 
kinetic energy is greater than said elemental breakpoint, calculating an unsealed flux wherein said 
unsealed flux equals A 3 ><EN YS1 ? wherein A 3 is a function of said elemental spectrum model and said 
baseline model, and wherein said A 3 is selected from the table: 



Baseline model 

Elemental 

spectrum model ^^^^^ 


worst day or peak flux 
in (cm 2 sr MeV/nuc)' 1 


worst week 

in (cm 2 sr MeV/nuc)' 1 


iron 


0.252948* 10 10 


0.2497 19xl0 9 


oxygen 


0.106702xl0 10 


0.667628*10 9 



and wherein y S i is a spectral index and is a function of said elemental spectrum model and said 
baseline model, and wherein said y S i is selected from the table: 



^^^^^ Baseline model 
Elemental ^ — 
spectrum model — 


worst day or peak flux 


worst week 


iron 


-4.52970 


-3.7610 


oxygen 


-4.14060 


-3.76850 



(I) in the case wherein the condition recited in step (G) is not satisfied, and wherein said 
kinetic energy is less than or equal to said elemental breakpoint, calculating an unsealed flux wherein 
said unsealed flux equals A 2 exp(-GxEn 1/4 )*En 1/4 , wherein A 2 is a function of said elemental 
spectrum model and said baseline model, and wherein said A 2 is selected from the table: 



Baseline model 



Elemental 
spectrum model 



worst day or peak flux 
in (cm 2 sr MeV/nuc)' 1 



worst week 

in (cm 2 sr MeV/nuc)"' 
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iron 


1.8991xl0 8 


3.0372x10 s 


oxygen 


4.9518xl0 8 


1.1307xl0 9 



and wherein G si is a spectral index and is a function of said elemental spectrum model and said 
baseline model, and wherein said G si is selected from the table: 



^ ^^^^ Baseline model 
Elemental ^^^^^^ 
spectrum model ^^^^^ 


worst day or peak flux 


worst week 


iron 


5.70 


5.70 


oxygen 


5.70 


5.70 



; and 



(J) in the case where said atomic number is between 3 and 92, inclusive, calculating a solar 
ion flux for said specified atomic number and said specified kinetic energy by multiplying said 
unsealed flux by a scale factor ratio, said scale factor ratio being the ratio of a scale factor for an 
element having said selected atomic number over a scale factor for said spectrum model element, 
wherein said scale factors are selected from the table: 
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ru 



Atomic 
Number 


Scale Factor 


Atomic 
Number 


Scale 
Faaor 


Atomic 
Number 


Scale 
Factor 


Atomic 
Number 


Scale 
Factor 


Atomic 
Number 


Scale 
Factor 


Atomic 
Number 


Scale 
Factor 


5 


0 


22 


4.377 
xlO" 3 


39 


4.878 
xlO* 6 


56 


4.878 
xlO" 6 


73 


2.195 
xlO" 8 


90 


4.878 
xlO" 8 


6 


4.704 
xlO' 1 


23 


4.088 


40 


1.22 
xlO" 5 


57 


4.878 
xlO" 7 


74 


2.439 
xlO" 7 


91 


0 


7 


1.2059 

xlO" 1 


24 


1.65 
xlO' 2 


41 


9.756 
xio 7 


58 


1.22 
xlO" 6 


75 


4.878 
xlO" 8 


92 


2.927 
xlO- 8 


8 


1 


25 


5.625 
xlO 3 


42 


4.878 
x 10" 6 


59 


1.951 
xlO 7 


76 


7.317 

xlO* 7 


3 


0 


9 


4.560976 

xlO 5 


26 


1 


43 


0 


60 


9.756 
xlO 7 


77 


7.317 

xlO* 7 


4 


0 


10 


2.1312 
xlO 1 


27 


1.303 
xio 2 


44 


2.195 

xlO" 6 


61 


0 


78 


1.463 
xlO* 6 






11 


1.744715 
x 1 0" 2 


28 


3.172 

xlO" 2 


45 


4.878 
xlO* 7 


62 


2.439 

xlO" 7 


79 


2.439 

xlO' 7 






12 


2.0624 

xlO' 1 


29 


3.048 

xlO" 4 


46 


.1.463 
xlO* 


63 


9.756 
xlO* 8 


80 


2.439 
xlO* 7 






13 


1.826829 


30 


7.457 

xlO" 4 


47 


4.878 
xlO* 7 


64 


4.878 
xlO' 7 


81 


2.195 
xlO 7 






14 


3.5935 
xio' 


31 


4.878 
xlO" 5 


48 


1.707 
xlO" 6 


65 


7.317 

xlO" 8 


82 


2.439 
xlO" 6 






15 


2.279675 
x 10" 4 


32 


1.22 
xlO" 4 


49 


2.195 

xlO* 7 


66 


4.878 

xlO" 7 


83 


1.463 
xlO 7 






16 


9.758 
xlO 2 


33 


7.317 

xlO* 


50 


4.878 
xlO* 6 


67 


9.756 
xlO* 8 


84 


0 






17 


1.680488 
xlO a 


34 


7.317 
xlO* 5 


51 


3.415 
xlO 7 


68 


2.439 
xlO 7 


85 


0 






18 


1.771545 

xlO" 3 


35 


9.756 
x I0" 6 


52 


7.317 
xlO" 6 


69 


4.878 
xlO 8 


86 


0 






19 


3.644715 
xlO* 4 


36 


4.878 
xlO* 5 


53 


1.463 
xlO" 6 


70 


1.951 

xlO' 7 


87 


0 






20 


4.826 
xlO* 2 


37 


7.317 
xlO* 6 


54 


6.585 
xlO" 6 


71 


4.878 
xlO" 8 


88 


0 






2! 


2.929 
xlO- 4 


38 


2.439 
xlO" 5 


55 


4.878 

xlO 7 


72 


1.951 
xlO 7 


89 


0 







wherein elements having atomic number of 20 or less are scaled to oxygen, and elements having 
higher atomic numbers are scaled to iron. 
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11. A method for determining a flux of solar ions of species having a range of atomic numbers 
between 3 and 92, wherein each of said species includes particles having a range of specified kinetic 
energies, for a satellite in a near earth orbit, comprising the steps: 

(A) for each kinetic energy of each of said species, performing steps (A) through (K) of claim 
10, and storing a flux value for each of said kinetic energies for each of said species; and 

(B) adding each of said flux values to obtain a total solar ion flux value. 

12. A method for modeling the effects of cosmic rays on microelectronics on a computer 
connected to the internet, comprising the steps: 

(A) receiving a login message from a remote user connected to the internet, and generating 
and transmitting back to said user with a script in response thereto an HTML main menu page, 
wherein said main menu page comprises prompts for said user to run one or more routines selected 
from the group consisting of: (1) a routine for calculating a geomagnetic transmission function for 
SEU-inducing particles, (2) a routine for calculating a flux of SEU-inducing particles in the near 
earth environment or in the environment shielded by earth's magnetosphere, (3) a routine for 
calculating a solid shielding transport function for SEU-inducing particles, (4) a routine for 
calculating a proton induced single event upset rate, (5) a routine for calculating a linear energy 
transfer rate, and (6) a routine for calculating a heavy ion induced single event upset rate, said main 
menu page further comprises prompts for said user to select a user request file for each of said 
routines; 

(B) in response to a message from said user to run at least one of said routines, using a user 
request file specified by said user, generating and executing with a script a system command for said 
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computer to run each of said specified routines; and 

(C) in response to the completion of all of said routines specified by said user, generating and 
transmitting to said user with a script an HTML page notifying said user of said completion. 

1 3 . The method of claim 1 2, wherein said main menu page generated and transmitted to said user 
further comprises prompts for said user to create or edit a user request file for any one of said 
routines, in lieu of running said one or more of said routines, and further comprising the step: 

(D) in response to a message from said user to create or edit a user request file, generating 
and transmitting to said user with a script an HTML page with prompts for the entry of fields for said 
user request file; and 

(E) in response to said user transmitting information for said user request file, storing said 
information in a user request file with a script. 

14. The method of claim 12, wherein said main menu page comprises prompts for said user to 
run a routine for calculating a geomagnetic transmission function for SEU-inducing particles, and 
wherein said step of generating and executing with a script a system command for said computer to 
run each of said specified routines comprises at least generating and executing a system command 
for said computer to run said routine for calculating a geomagnetic transmission function for SEU- 
inducing particles, in accordance with claim 1 

15. The method ofxlaimJ 2, wherein said main menu page comprises prompts for said user to 
run a routine for calculating a flux of SEU-inducing particles in the near earth environment or in the 
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environment shielded by earth's magnetosphere, and wherein said step of generating and executing 
with a script a system command for said computer to run each of said specified routines comprises 
at least generating and executing a system command for said computer to run said routine for 
calculating said flux of SEU-inducing particles in the near earth environment or in the environment 
shielded by earth's magnetosphere, wherein said routine comprises calculating a flux of solar heavy 
ions in accordance with claim 10. 
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